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ABSTRACT 

The large majority of EGRET point sources remain without an identified low-energy counterpart, 
and a large fraction of these sources are most likely extragalactic. Whatever the nature of the extra- 
galactic EGRET unidentified sources, faint unresolved objects of the same class must have a contribu- 
tion to the diffuse extragalactic gamma-ray background (EGRB). Understanding this component of 
the EGRB, along with other guaranteed contributions from known sources, is essential if we are to use 
this emission to constrain exotic high-energy physics. Here, we follow an empirical approach to esti- 
mate whether a potential contribution of unidentified sources to the EGRB is likely to be important, 
and we find that it is. Additionally, we show how upcoming GLAST observations of EGRET unidenti- 
fied sources, as well as of their fainter counterparts, can be combined with GLAST observations of the 
Galactic and extragalactic diffuse backgrounds to shed light on the nature of the EGRET unidentified 
sources even without any positional association of such sources with low-energy counterparts. 
Subject headings: gamma rays: observations - gamma rays: theory - diffuse radiation 



1. INTRODUCTION 

The origin of the isotropic diffuse emission (Sreekumar 
et al. 1998) in energies between 100 MeV and 20GeV, 
detected by the Energetic Gamma-Ray Experiment Tele- 
scope (EGRET) aboard the Compton Gamma-Ray Ob- 
servatory, remains one of the great unknowns of GeV- 
energy astrophysics. There are two major questions that 
still remain unanswered. 1. How much of the diffuse emis- 
sion detected by EGRET is, in fact, extragalactic, and 
what is the spectrum of this extragalactic background? 
And 2. what fraction of the extragalactic emission can 
be attributed to each of the observationally established 
classes of gamma-ray emitters? Despite the associated 
uncertainties, these two issues are critical in any attempt 
to use gamma-ray observations to constrain exotic high- 
energy physics and yet-undetected classes of theorized 
gamma-ray emitters. 

To answer the first question, a good understanding of 
the Galactic diffuse emission in the EGRET energy range 
is required. In order to obtain the intensity and spec- 
trum of the extragalactic emission from the EGRET sky 
maps, the Galactic emission needs to be modeled and 
subtracted. This is made complicated by the discrep- 
ancy between the observed Milky Way spectrum in en- 
ergies of > lGeV and theoretical expectations (Hunter 
et al. 1997). The observed spectrum is more shallow 
than model predictions based on the local demodulated 
cosmic ray spectrum. This deviation is known as the 
"GeV excess", and although various explanations have 
been proposed to account for part or all of the discrep- 
ancy, its origin remains a matter of debate (e.g., Pohl et 
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al. 1997; Busching et al. 2001; Strong et al. 2004b; Ka- 
mae et al. 2005; de Boer et al. 2006; Strong 2006; Stecker 
et al. 2007). As a result, determinations of the gamma- 
ray background using different Galactic emission models 
yield very different answers, both in intensity and in spec- 
trum, despite being based on the same set of observations 
(e.g., Sreekumar et al. 1998; Strong et al. 2004a). 

Attempts to answer the second question have been 
plagued by uncertainties in the cosmic density and evo- 
lution of the two established classes of extragalactic 
gamma-ray emitters: normal galaxies and blazars. Our 
observational knowledge of the gamma-ray properties of 
normal galaxies is very limited, as the sample of normal 
galaxies which have been observed in gamma rays con- 
sists of only the Milky Way and a marginal detection 
of the Large Magellanic Cloud (Sreekumar et al. 1992; 
Hunter et al. 1997; Hartman et al. 1999). For this reason, 
the accuracy of theoretical estimates of the contribution 
of normal galaxies to the gamma-ray background is un- 
avoidably at the order-of-magnitude level (e.g., Lichti et 
al. 1978; Pavlidou & Fields 2002). But even in the case 
of blazars, which are by far the most numerous and best 
studied class of identified gamma-ray emitters, estimates 
of their contribution to the gamma-ray background vary 
from a few percent to 100% of the background origi- 
nally reported by the EGRET team (e.g., Padovani et 
al. 1993; Stecker & Salamon 1996a; Kazanas & Perlman 
1997; Mukherjee & Chiang 1999; Miicke & Pohl 2000; 
Narumoto & Totani 2006; Dermer 2007). 

The issue is further complicated by the existence of 
171 sources which, at the time of publication of the 3rd 
EGRET catalog (hereafter 3EG; Hartman et al. 1999), 
had not been positively or potentially associated with 
a lower-energy counterpart. These sources are collec- 
tively known as the unidentified EGRET sources, and 
they are more numerous than any established group of 
gamma-ray emitters. The distribution of these sources 
on the sky is such that a Galactic feature can be clearly 
distinguished - however a large number of sources are 
located away from the Galactic plane and the Galac- 
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tic center 6 . No more than a handful of sources can be 
associated with the Milky Way halo if the Milky Way 
is not many times brighter in gamma-rays than simi- 
lar galaxies such as M31 (Siegal-Gaskins et al. 2007). 
Hence, it is almost certain that the EGRET unidenti- 
fied sources include a significant extragalactic compo- 
nent. Although the nature of these extragalactic sources 
remains unknown, it is reasonable to believe that there is 
a large number of fainter, unresolved objects of the same 
class, which are guaranteed to have some contribution 
to the extragalactic gamma-ray background (EGRB). If 
these sources represent yet unidentified members of some 
known class of gamma-ray emitters (e.g. blazars), then 
excluding them from any calculation of the contribution 
of the parent class to the diffuse background would lead 
to a significantly underestimated result due to an incor- 
rect normalization of the bright-end of the gamma-ray 
luminosity function. If they represent an unknown class 
of gamma-ray emitters, then the contribution of their un- 
resolved counterparts to the diffuse emission would sig- 
nificantly limit the diffuse flux left to be attributed to 
known classes, exotic processes, and truly diffuse emis- 
sion. 

Hence, some contribution of unresolved unidentified 
sources to the EGRB is certain. It is therefore clear 
that until we either answer the question of the nature of 
unidentified sources or derive some strong constraint in- 
dicating that a possible contribution of such unresolved 
objects to the EGRB would indeed be minor, we cannot 
hope to fully understand the origin of the EGRB. 

Detailed predictions for the level of the unidentified 
source contribution to the EGRB involve important un- 
certainties: since no low-energy counterparts have been 
identified, we have no estimates of distance, and there- 
fore no estimates of the gamma-ray luminosities of these 
sources. As a result, very few constraints can be placed 
on their cosmic distribution and evolution. However, 
very simple estimates can offer some guidance on whether 
ignoring this EGRB component may be a safe assump- 
tion to make. 

For example, we can use the number of unidenti- 
fied sources, the minimum flux resolvable by EGRET, 
and the observed intensity to the extragalactic gamma- 
ray background to place rough limits on the distance 
scales associated with resolved and unresolved uniden- 
tified sources so that unresolved sources do not overpro- 
duce the background. A population of unbeamed, non- 
evolving, single-luminosity sources uniformly distributed 
in Euclidian space are resolvable out to a distance D 
by an instrument of number flux sensitivity F m i n . The 
relation between D, F m i n , and number luminosity L in 
this case is simply L — 4nD 2 F m i n . If the instrument de- 
tects N such sources, their number density n SOU rcc can 
be estimated to be n sourco = 3N/4TrD 3 . If the same dis- 
tribution of sources continues out to a distance d > D, 
the isotropic intensity (photons per unit area per unit 
time per unit solid angle) from the unresolved mem- 
bers of this population will be I U nrcs = f D dl s hcii, where 
dlshc\\ = (l/47r)(n sourcc 47rr 2 (ir)L/(47rr 2 ) is the contribu- 
tion from sources within a spherical shell located at a 

6 Note however that the presence of sources at high latitudes does 
not, in itself, constitute proof that these objects are extragalactic 
(see e.g. the Gould Belt discussion in Gehrels et al. 2000). 



distance r from the observer. Substituting our results 
for n source and L above, and performing the integral, we 
obtain 

/unrcs = 3NF min (d - D)/4nD . (1) 

If we require that the unresolved emission from this pop- 
ulation does not exceed the EGRB observed by EGRET 
(^unros < /egrb), we obtain Zegrb > 3NF min (d - 
D)/4itD. Substituting N ~ 100 for the population of 
extragalactic unidentified sources (see discussion in 
-Fmin ~ 10- 7 phcm 2 s- 1 for the sensitivity of EGRET, 
and Iegrb ~ 10 -5 phem 2 s _1 sr^ 1 (Sreekumar et al. 
1998), we obtain d < 6D. This result implies that the 
largest distance out to which such a distribution of ob- 
jects persists cannot be larger than a few times the dis- 
tance out to which these objects are currently resolved, 
since in any other case these sources would overproduce 
the EGRB. It is therefore conceivable that the uniden- 
tified source contribution to the EGRB is significant, if 
not dominant. 

In this work, we approach the problem from a purely 
empirical point of view. Instead of attempting to predict 
the level of a diffuse component due to unresolved objects 
of the same class as unidentified EGRET sources, we try 
to assess whether there are any empirical indications that 
this component is, in fact, minor. We construct samples 
of unidentified sources which, based on their sky distribu- 
tion, are likely to consist mostly of extragalactic objects. 
Under the assumption that the majority of these sources 
can be treated as members of a single class of gamma-ray 
emitters, we seek to answer the following three questions: 

(1) Is it likely that unresolved objects of the same class 
could have a significant contribution to the EGRB at 
least in some energy range? 

(2) How would the collective spectrum of their emis- 
sion compare to the measured spectrum of the EGRB 
deduced from EGRET observations? 

(3) How are GLAST observations expected to improve 
our understanding of the nature of unidentified sources, 
based on the insight gained from our analysis? 

This paper is structured as follows. In fJ3] we discuss 
the samples of resolved unidentified sources used in our 
analysis. Our formalism for constructing the collective 
emission spectrum of unresolved unidentified sources is 
presented in $3] Inputs from EGRET data used in our 
analysis are described in In Sj5] we describe our re- 
sults, and in $6]we discuss how these results are expected 
to improve once GLAST observations become available. 
Finally, we summarize our conclusions in <J7] 

2. SOURCE SAMPLES 

The observational input constraining our calculations 
will be the fluxes of the unidentified sources and the spec- 
tral index of each source. The fluxes we use for our anal- 
ysis are the P1234 photon fluxes quoted in 3EG. These 
fluxes were determined from the combined maps from all 
viewing periods during which any particular object was 
within the EGRET field of view. The spectral index a 
of each object is also taken from 3EG, and was derived 
through a fit assuming that the photon flux of each ob- 
ject obeys F1234 = Fq(E / Eo)~ a where E is the photon 
energy. 

We seek to constrain a possible contribution to the 
EGRB from objects of the same class as unidentified 



3 



90 




''■•-•'"■&;;;^v,.j.. 1 v:i;;V^^ v "' 

-90 



Fig. 1. — Source exclusion mask that has been applied in our 
analysis. Sources with coordinates within the shaded regions have 
been excluded from our analysis. The locations of sources in our 
samples are also indicated. X: possible 3EG blazars (included 
in Sample 1 only); □: originally unidentified in 3EG, with later 
suggested associations (included in Samples 1 and 2); A: sources 
included in all three samples (unidentified in 3EG, no counterpart 
suggested to this day). 



sources. For this reason, we would like to base our cal- 
culations on a resolved source sample which is as repre- 
sentative as possible of the population of extragalactic 
unidentified EGRET sources. However, the population 
of sources which, at the time of publication of 3EG, did 
not have confident associations with known astrophysical 
objects, clearly exhibits a Galactic feature, with many of 
the sources concentrated around the Galactic plane and 
at small Galactic longitudes (close to the Galactic cen- 
ter). Thus, the original population of EGRET unidenti- 
fied sources is not an appropriate sample as it includes 
many sources which are most likely Galactic. 

Since these sources are unidentified, there is no coun- 
terpart information which could be used to distinguish 
between Galactic and extragalactic objects. Addition- 
ally, inferences based on their spectral indices and vari- 
ability are also very hard to make due to the large error 
bars associated with measurements of these quantities 
(see Siegal-Gaskins et al. 2007). Finally, both Galactic 
and extragalactic identified objects exhibit a large range 
of luminosities, so that inferences about the location of 
unidentified sources based on their fluxes are not possible 
either. The only relatively robust indication concerning 
the location of these sources is their position on the sky. 
A Galactic population associated with the disk and bulge 
will most likely be concentrated around low Galactic lati- 
tudes and small Galactic longitudes. On the other hand, 
a halo population, which may appear more uniformly 
distributed on the sky, cannot have too many members 
among the resolved EGRET unidentified sources, based 
on the gamma-ray properties of galaxies similar to the 
Milky Way (Siegal-Gaskins et al. 2007). Thus, excluding 
objects around the Galactic plane and the Galactic cen- 
ter will most likely minimize any Galactic contamination 
of the sample of unidentified sources. 

For this reason, we have excluded from our studies 
all sources with Galactic latitude |6| < 10° for Galac- 
tic longitudes \l\ > 40° and Galactic latitude |b| < 30° 
for Galactic longitudes |/| < 40°. The source exclusion 
map we have applied is plotted in Fig. [TJ and is identical 



to the mask applied by Sreekumar et al. (1998) in their 
selection of EGRET diffuse emission data to use for the 
derivation of the EGRB. The solid angle associated with 
this mask is r2 mas k ~ 7r, leaving 37r of allowed solid 
angle. For this reason, and under the assumption that 
an extragalactic population will be uniformly distributed 
on the sky, we have renormalized all derived flux distri- 
butions by a factor of 4/3. 

Additionally, we have excluded from all samples ob- 
jects J0516+2320, J0852-1216, J1424+3734, and J2241- 
6736 because no P1234 flux was quoted for them in 3EG, 
and objects J0824-4610, J0827-4247, J0828-4954, J0841- 
4356, J0848-4429, and J0859-4257 because they were 
marked as possible or likely artifacts in 3EG. 

We have applied further selection criteria which result 
in three distinct samples of resolved unidentified sources. 
The use of different samples will allow us to test the 
sensitivity of our results to the details of the included 
members. 

Sample 1 includes all sources which were unidentified 
at the time 3EG was published and all sources which 
were characterized as "possible blazars" (identification 
code "a" in 3EG). These latter objects are routinely not 
included in studies of the blazar luminosity function, so 
it is reasonable to examine the potential contribution of 
this class to the EGRB by including it in our analysis. 
Sample 1 includes 84 sources. 

Sample 2 only includes sources which were unidenti- 
fied at the time 3EG was published. In this case, we do 
not consider "possible blazars" to belong to this sample. 
Sample 2 includes 63 sources, and is a subset of Sample 
1. 

Sample 3 is the most conservative of our samples, in 
that it contains the smallest number of objects. In this 
case, not only have we excluded all "probable blazars" , 
but also any object for which even a tentative or potential 
identification has been suggested since the publication of 
3EG 7 . Sample 3 includes 53 sources, and is a subset of 
Sample 2. Sources which are present in Sample 2 but not 
in sample 3 are the following: J0010+7309, J0724-4713, 
J1212+2304, J1825+2854, J1835+5918, J0329+2149, 
J0725-5140, J1249-8330, J1621+8203, and J1959+6342. 
The locations of all the sources in our samples are shown 
in Fig. CD 

3. FORMALISM 

A quantity which, under the assumption of isotropy, 
fully describes the EGRB is the differential photon in- 
tensity Ie (photons per unit area-time-energy-solid an- 
gle). In this section, we calculate Ie for a population of 
unresolved sources. We make the following assumptions 
for the population of unresolved gamma-ray sources we 
consider. 

1. Each of the sources has an energy spectrum which, 
in the EGRET and GLAST energy ranges, can be 
well described as a single power law of spectral 
index a, and a total photon flux F& at energies 

7 Sample compiled and actively maintained by C. Brown and 
available online at http://GeVsky.org We stress that suggested 
associations with lower-energy counterparts included in this com- 
pilation are simply results reported in recent publications; no effort 
was made to evaluate the significance and validity of these coun- 
terparts by a single uniform standard. 
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E > 100 MeV (where Fg is the photon flux in units 
of 10~ 8 photons cm~ 2 s _1 ). 

The cumulative flux distribution, N(> Fg) (num- 
ber of sources with flux greater than Fg versus F 8 ), 
of this population can be adequately described, be- 
tween fluxes F 8m j n and F g . max , as a power law, 



N(> F 8 ) = CF 8 



(2) 



3. The probability distribution of spectral indices of 
the population p(a) is independent of source lumi- 
nosity and does not evolve with redshift. 

The differential flux of each object of this population, 
Fe, is related to F 8 through 
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provided a > 1. Each unresolved source with a flux Fg 
at energies > Eq, has a contribution Ie,i{F»,ol) to the 
diffuse emission which is given by 
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where a is the spectral index of the source, and Eq = 
100 MeV in the case of 3EG fluxes. The 1 /An normaliza- 
tion factor comes from assuming an isotropic distribution 
of sources, the collective emission of which is uniformly 
distributed over the celestial sphere. 
The differential flux distribution of the population is 

'' - V kCF 8 -^ . (5) 



dF 8 

The collective diffuse emission due to this population of 
unresolved sources with fluxes between Fg im ; n and -Fs.max 
will be: 
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where I = 10 _7 cm -2 s -1 sr _1 GeV _1 /(47i")- Physically, 
the upper-limit flux i^.max represents the sensitivity 
limit of the telescope (objects with flux higher than this 
are resolved and do not contribute to the background), 
while the lower-limit flux F 8:m i n represents the flux be- 
low which the approximation of Eq. ([2|) for the cumula- 
tive flux distribution breaks down. The features of the 
spectral shape of this emission in the simple case when 
p{a) is a Gaussian are discussed in Appendix [Ul 

4. INPUTS FROM EGRET DATA 

Equation ^ can provide a first simple estimate for the 
possible contribution of unresolved sources of the same 
class as extragalactic unidentified EGRET sources to the 
isotropic diffuse gamma-ray background. The required 
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Fig. 2. — Cumulative flux distribution of the resolved unidentified 
sources. Solid line: sample 1; dotted line: sample 2; dashed line: 
sample 3. The thin solid lines represent the location of the median 
value of sample 1 for fluxes below 6 X 10 — 8 photons cm~ 2 s~ 1 . The 
vertical lines correspond to the lowest flux to which we can extrapo- 
late a power-law fit to the cumulative flux function of each sample 
in the flux range Fg max < Fg < Fg jv = s, before the cumulative 
emission from unresolved point sources exceeds the measurement 
of the EGRB at low energies. The light grey vertical line shows the 
lowest flux to which we can perform the same extrapolation if the 
flux function of sample 1 combined with identified blazars outside 
our mask is used instead. 



inputs are the cumulative flux function of Eq. @ , and the 
spectral index distribution of the source class, p(a). In 
this section, we derive these inputs from the three source 
samples discussed in fJ5] Possible source variability has 
not been accounted for in this analysis. A discussion 
of possible effects of source variability, as well as our 
reasoning for not considering them for the purposes of 
our study, are presented in Appendix [Bl 

4.1. The cumulative flux distribution 

We use data from the resolved objects in our three 
samples to construct the cumulative flux distribution of 
EGRET unidentified sources. The cumulative flux dis- 
tributions of the three samples are shown in Fig. (J5]). 
Due to non-uniformities in EGRET exposure, the lowest 
resolvable flux is not a hard-set number, but rather the 
efficiency in resolving low-flux objects decreases gradu- 
ally as the flux decreases, which is one of the major fac- 
tors in the slow turnover of the cumulative distribution 
function. For this reason, we only use data within the 
range Fg imax < Fg < Fgjv=5 to obtain a fit to the cumu- 
lative flux distribution. We take Fg jmax = 6, and assume 
that all sources with fluxes above this limit have been 
resolved by EGRET. We define Fg^^s to be the flux for 
which, in each sample, 5 objects with fluxes higher than 
or equal to it have been resolved. The best-fit values of 
C and k for each sample are shown in Table [T] As we 
can see in Fig. © as well as from the parameter values in 
Table [TJ the slopes of the cumulative flux distribution in 
the regime where most objects are expected to have been 
resolved by EGRET are fully consistent with each other. 
This in turn implies that the data are consistent with the 



hypothesis that most members of all samples have been 
drawn from a single population of extragalactic emitters. 

However, we note that the flux distribution slope, 
close to ~2.5 in all cases, is steeper than the slope that 
would be expected from a single-luminosity population 
of sources with a uniform distribution in a flat cosmol- 
ogy (which is equal to 1.5, e.g. Dermer 2007), and is 
also steeper than the slope of the flux distribution of 
the confidently identified blazar population. This result 
may reflect the cosmological distribution and evolution 
properties of the population of extragalactic unidentified 
sources. However, it may also originate in a selection 
effect in the identification of gamma-ray sources: the po- 
sitional identification of brighter sources is more frequent 
as more photons are detected from these sources which 
allows for a more accurate pinpointing of their location. 
In addition, multiwavelength campaigns for the identifi- 
cation of gamma-ray sources naturally target the bright- 
est objects first. For these reasons, the high-flux end 
of the flux distribution may be preferentially depleted, 
which would lead to an apparent steepening of the dis- 
tribution of the sources that remain unidentified. Finally, 
the steepness of the slope could be an effect of viewing a 
distribution with real curvature in a very small dynam- 
ical range in flux. A simple example quantifying such 
a possible effect is the following: if we were to treat all 
sources in Sample 1 and all confidently identified blazars 
which are located outside our exclusion mask as a single 
population, then the cumulative flux distribution power- 
law fit in the flux range 6 < F$ < l<8,;v=5 would have a 
slope k — 1.63 ± 0.03, much closer to 1.5. 

Finally, we should add a cautionary note on fits to 
the cumulative, rather than differential, flux distribution 
function, such as the one presented here. Our fits purpos- 
fully do not account for uncertainties in each bin since, 
from a statistical point of view, such uncertainties are 
almost perfectly correlated (each bin contains the same 
data as its adjacent ones plus/minus one data point). On 
the other hand, the poor dynamical range in flux and 
the dependence of uncertainties of the differential flux 
function on the size of the bin make constructing and 
fitting the differential flux function problematic as well. 
For these reasons, it is important to not over-interpret 
these results, which should be viewed as broad, order-of- 
magnitude assessments of the behavior of the flux dis- 
tribution, rather than robust statistical evaluations and 
strict constraints. 

4.2. The spectral index distribution 

The second important observational input in our cal- 
culation is the distribution of spectral indices of uniden- 
tified objects p{a). In the limit that the spectral index 
of sources is independent of source redshift and luminos- 
ity, the spectral index distribution uniquely determines 
the spectral shape of the collective emission from the 
unresolved population. The spectral shape of the unre- 
solved emission provides, in turn, an additional tool to 
assess the possibility that unidentified sources constitute 
a dominant contribution to the EGRB through shape 
comparison between the observed and the predicted spec- 
tra of diffuse extragalactic emission. 

In our analysis we adopt the assumption that the spec- 
tral index distribution of unresolved unidentified sources 
is the same as that of the resolved unidentified sources. 
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Fig. 3. — Dashed line: Sreekumar et al. (1998) determination of 
the EGRB. Solid line: Strong et al. (2004a) determination of the 
EGRB (best guess). Dotted lines: Strong et al. (2004a) determi- 
nation of the EGRB (systcmatics- based limits). Thick solid line: 
collective spectrum shape of unresolved unidentified sources us- 
ing maximum-likelihood intrinsic spectral index distribution (ISID) 
parameters for Sample 1 (this work). The amplitude of the spec- 
trum does not constitute a prediction; rather, it has been maxi- 
mized so that it does not exceed the observed EGRB in any energy 
bin above 100 MeV. Thick grey lines: lc uncertainty in the spec- 
tral shape of unresolved unidentified source emission for Sample 1, 
entering through the statistical uncertainty in the ISID parameters. 



The latter can be deduced from measurements of the 
spectral index a for sources in each one of our samples, 
following the method presented in Venters & Pavlidou 
(2007). The details of this calculation, properly account- 
ing for measurement uncertainties in individual spectral 
indices of resolved sources, are presented in Appendix lAl 

5. RESULTS 

As is immediately obvious from Eq. © , the shape and 
the overall normalization Iq of the cumulative emission 
from unresolved unidentified sources are decoupled under 
our assumptions. The shape of the spectrum depends on 
the spectral index distribution, while the normalization, 
given a fit to the flux distribution, depends only on i^min 
(the value of the flux where the extrapolated power law 
breaks). 

Assuming that, close to the EGRET flux limit, the 
flux distribution does not evolve drastically, an extrapo- 
lation of the measured flux function to lower fluxes can be 
considered representative of its behavior in the low-flux 
regime. This assumption is less likely to hold as the lim- 
iting flux to which we are extrapolating becomes lower: 
the flux distribution will eventually exhibit a break due 
to cosmological effects and/or luminosity evolution. The 
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question is how far in the low-flux regime our extrapo- 
lation must continue before we get a significant contri- 
bution of unresolved unidentified sources to the gamma- 
ray background. If the answer is "not very far", then 
our extrapolation may indeed be reasonable over such 
a small range, and unresolved, unidentified sources are 
likely make up a considerable fraction of the EGRB. On 
the other hand, if we need to extrapolate the flux distri- 
bution down to fluxes very low compared to the resolved 
flux range, then it is quite unlikely that our extrapola- 
tion is valid throughout the flux regime in which we use 
it. In such a case it is doubtful that the actual flux dis- 
tribution of the unresolved unidentified sources is such 
that unidentified sources make up a significant portion 
of the EGRB. 

We use observations of the EGRB to determine by 
what value of -Fs.min the extrapolated flux distribution 
power-law must break, so that the background is not 
overpredicted at any energy range. To ensure that we 
are placing the most conservative constraint (i.e. that 
we allow extrapolation to the lowest fluxes possible), we 
will use the most generous observational determination 
of the EGRB, which is that of Sreekumar et al. (1998). 
The requirement that this background is not exceeded for 
any energy above 100 MeV returns the limiting breaking 
fluxes indicated by the three vertical lines in Fig. [2] Our 
results then indicate that, for any of our three samples, 
we would only need to extrapolate the cumulative flux 
distribution for at most an order of magnitude below 
the lower limit of the resolved flux range to have unre- 
solved unidentified sources comprise most of the EGRB, 
at least at low energies. This is not an extreme extrap- 
olation, and therefore a significant contribution by the 
"unidentified" class to the diffuse background is likely. 

In order to interpret this result correctly, we should 
assess the qualitative and quantitative uncertainties as- 
sociated with it. First of all, a power-law extrapolation 
of the flux distribution to lower, unresolved flux values 
even for only one order of magnitude in flux is by no 
means self-evident, although it is the simplest and most 
straight-forward assumption. For example, the Dermer 
(2007) fits to the flux distribution of resolved blazars 
(using a model based on the cosmological evolution of 
black hole jet sources) increase less steeply than a power 
law even very close to the EGRET sensitivity limit. It is 
conceivable that unresolved members of the extragalactic 
unidentified source class could have a comparable cosmic 
evolution, and a similar behavior in the flux function. 
In fact, our tight constraints on the power-law extrap- 
olation of the flux function may be taken to imply just 
that: the cosmological evolution of these sources has to 
be such that the flux function deviates from a power-law 
form very fast, to ensure than unresolved members of 
this class do not overproduce the EGRB. 

Second, statistical and systematic uncertainties in the 
power-law fits to the cumulative flux function have to 
be taken into account. In our case, the systematic un- 
certainties associated with the sample selection and the 
small dynamical range in flux over which the flux func- 
tion is sampled dominate over the statistical uncertain- 
ties in the power-law slope derived in the selected flux 
range, shown in Table [TJ The most extreme case is the 
one where unidentified sources are considered as a unified 
sample with identified blazars outside the mask. In this 



case, as discussed in £|4. 1 1 the flux function slope is con- 
siderably shallower, allowing for a more extended extrap- 
olation to low fluxes before the EGRB is overproduced. 
The lowest flux limit allowed for the extrapolation in this 
case is shown with the grey line in Fig. [21 and is a moder- 
ate factor ~ 3 below the limits of the unidentified source 
samples when considered alone. 

A comparison between the shape of the cumulative 
emission spectrum of unresolved unidentified sources and 
the EGRET EGRB is shown in Fig. H The dashed line 
shows a single power-law fit to the Sreekumar et al. 
(1998) determination of the EGRB. The thin solid line is 
the more recent redetermination of the gamma-ray back- 
ground by Strong et al. (2004a), in which they used a 
more detailed model of the Milky Way diffuse emission 
to subtract the Galactic component from the EGRET 
diffuse sky map. The dotted lines are the systematic un- 
certainties in the EGRB determination of Strong et al. 
(2004a), entering through their model of the Galaxy. Our 
calculation of the spectrum of the unresolved unidenti- 
fied source component is shown with the thick solid line. 
The maximum-likelihood parameters from Sample 1 were 
used for p(a) i n th is calculation, although, as becomes 
clear from Fig. IA2I in the Appendix, the spectral shapes 
resulting from all three samples are consistent with each 
other. At low energies, where the systematics are low, 
the spectrum of the unidentified component is in excel- 
lent agreement with the EGRB observational spectrum 
of Strong et al. (2004a). At higher energies, where the 
systematics are large, the unidentified component spec- 
trum is largely within systematics except at very high 
energies. If unidentified sources are indeed a dominant 
contribution at relatively low energies, then this result 
may be perceived as a tantalizing hint that at the highest 
energies of the EGRET range a new type of contribution, 
(e.g., from high-energy peaked BL Lacs, or from annihi- 
lating dark matter) may become important at ~ 1 GeV. 

It should be noted that the amplitude of the cumu- 
lative spectrum plotted in Fig. [3] does not constitute a 
prediction; rather, it is the maximum possible emission 
allowed from unresolved unidentified sources so that the 
observed EGRB is not exceeded in any energy bin above 
100 MeV. In Fig. [3]we also overplot, with the thick grey 
lines, the la uncertainty in the spectral shape of the 
unresolved emission from unidentified sources, entering 
through uncertainties in the ISID parameters. 

6. PROSPECTS FOR GLAST 

The launch of GLAST in 2008 will provide significant 
new insight into the nature of unidentified sources and 
their possible contribution to the EGRB. The ideal so- 
lution to the unidentified source puzzle would be, of 
course, the direct positional association of all uniden- 
tified sources with undisputed low-energy counterparts. 
This would then allow us to build more confident mod- 
els for the unresolved members of these classes of ob- 
jects. However, such an outcome is unlikely, as the large 
number of possible counterparts and the large number of 
sources which we expect GLAST will be able to resolve 
make multi-wavelength campaigns for every single source 
impractical. 

However, there is another definitive test that GLAST 
will be able to perform which does not require confi- 
dent identification of each source to provide information 



7 



about the likely nature of unidentified sources as a pop- 
ulation. With the increased flux sensitivity of GLAST, 
many more objects of the same class will be resolved. If 
these objects are mostly extragalactic, as assumed here, 
and if they comprise a significant fraction of the EGRET 
EGRB, as suggested by our analysis, then there will be 
an associated decrease of the EGRB from its EGRET lev- 
els, equal to the all-sky-averaged intensity of the newly 
resolved objects. The flux sensitivity of GLAST is ex- 
pected to be about 50 times better than that of EGRET 8 . 
Therefore GLAST will be able to probe the flux distribu- 
tion of unidentified sources down to fluxes close to i^mm 
and definitively test our empirical estimate. 

If, on the other hand, even high-latitude unidentified 
sources are mostly Galactic (which would imply that the 
Milky Way may be unusually bright in gamma-rays), 
then there will be an associated reduction of the Milky 
Way diffuse emission rather than of the isotropic back- 
ground (for a discussion on a possible contribution of 
a large number of Galactic point sources to the Galactic 
diffuse emission and the role of such sources in explaining 
at least in part the origin of the GeV excess, see Strong 
2006). 

GLAST will not only represent a leap in instrument 
sensitivity, but also in energy range, as it will allow access 
to the yet-unexplored range between a few tens and a few 
hundreds of GeV. Sources with fluxes high enough that 
measurements of their spectra can be extended in this 
range will have an additional indicator of their location 
(Galactic vs extragalactic): if they are located at cosmic 
distances, a feature generated by attenuation of gamma- 
rays on the extragalactic background light may be de- 
tectable; on the other hand, if they are located within 
the Milky Way, no such feature should be present. 

In addition, GLAST observations will improve the in- 
puts used in our calculation. Spectral indices, at least for 
the brighter sources, will be measured with much greater 
accuracy. This will improve our constraints on p(a) and 
the shape of the cumulative emission from unidentified 
sources. Additionally, spectral discrimination between 
different classes of unidentified emitters may become pos- 
sible. Finally, the dynamical range over which fluxes will 
be measured will increase significantly, and any biases 
in the calculation of the slope of the flux distribution of 
high-latitude unidentified sources will become less pro- 
nounced. 

7. DISCUSSION AND CONCLUSIONS 

In this work, we have used a purely empirical model 
to explore the possibility that unresolved gamma-ray 
sources of the same class as unidentified EGRET sources 
have an appreciable contribution to the EGRB. We have 
argued that some unidentified source contribution to the 
gamma-ray background is guaranteed. We have addi- 
tionally found that (1) if most high- latitude unidentified 
sources are assumed to be extragalactic, a one order of 
magnitude extension of the cumulative flux distribution 
to lower energies without breaks implies a significant con- 
tribution to the EGRB, at least at the lower part of the 
EGRET energy range; and (2) the spectrum of the cu- 
mulative emission of such unresolved sources would be 
very consistent with the observational determination of 

8 http:/ /www-glast. stanford.edu/mission. html 



Strong et al. (2004a) of the EGRET EGRB within sys- 
tematics. 

We emphasize that the purpose of this study is not to 
estimate, even at the order-of-magintude level, the dif- 
fuse flux expected from extragalactic unresolved uniden- 
tified sources, but rather to place constraints on the flux 
distribution of these objects under the constraint that 
the measured background should not be exceeded in any 
energy interval. Our treatment is therefore different in 
purpose and spirit from most past work aiming to es- 
timate the level of the contribution of different popula- 
tions to the extragalactic gamma-ray background. Such 
work has been traditionally of two types. The first in- 
volves population models built from some understanding 
of the physics of the sources (such as, e.g., the models 
of Miicke & Pohl 2000 and Dermer 2007 in the case of 
blazars; Miniati 2003, Keshet et al. 2003, Blasi, Gabici, & 
Brunetti 2007 in the case of clusters of galaxies; Lichti, 
Bignami & Paul 1978 and Pavlidou & Fields 2002 in 
the case of normal galaxies) . The second involves mod- 
els of the population luminosity function based on our 
knowledge of the source population from other wave- 
lengths and normalized to fit EGRET data (such studies 
require a sample of detected, identified members, and 
are therefore applicable to blazars only, e.g. Chiang et 
al. 1995, Stecker & Salamon 1996, Mukherjee & Chiang 
1999, Narumoto & Totani 2006). In our case, lacking 
any knowledge of the physics of sources as well as of 
even the bright end of the luminosity function (since in 
absence of identifications no redshift and hence no lumi- 
nosity can be derived for any of the sources), we have re- 
versed the problem. Instead of using some assumed cos- 
mic evolution for the sources to derive the expected level 
of contribution to the EGRB, as in all of the investiga- 
tions mentioned above, we have used the tightest possible 
constraint on the allowable EGRB contribution of these 
sources (the observed EGRET background) to constrain 
the flux distribution of the unresolved sources. Our con- 
clusions for the overall expected unresolved unidentified 
source intensity come from the observation that our con- 
straints on the flux distribution are indeed very tight; 
hence, the contribution of the unidentified sources to the 
EGRB is likely to be high. 

Our analysis suggests that any model of the EGRB 
would be incomplete without some treatment of the 
unidentified source contribution. The results of our em- 
pirical model therefore motivate the pursuit of specific 
population models for the unidentified sources. Although 
such models involve a more restrictive set of assump- 
tions and increased uncertainties, they can provide more 
concrete predictions for the luminosity function of unre- 
solved objects. Once a luminosity function model is as- 
sumed, and under the assumption that unresolved mem- 
bers of this class do indeed contribute most of the ex- 
tragalactic diffuse emission, additional estimates can be 
made regarding the distance scales associated with this 
population. The dynamical range in fluxes of the re- 
solved members of the population, as can be seen in Fig. 
[21 is smaller than an order of magnitude. For a single- 
luminosity population, this corresponds to a factor of 3 
dynamical range in distance, which increases or decreases 
if the typical source luminosity increases or decreases 
with increasing redshift respectively. Such arguments 
can constrain not only the distance scales, but also the 



number density, intrinsic brightness, and evolution of the 
resolved and unresolved objects, if in a model-dependent 
fashion. Additional constraints and predictions, once a 
luminosity function has been assumed, can be derived 
through the multi-messenger and multi-wavelength ap- 
proach. If the emission from extragalactic unidentified 
sources is primarily of hadronic origin, it will be ac- 
companied by neutrinos at comparable fluxes, and that 
may have potentially observable consequences in the TeV 
range for future km 3 neutrino detectors such as IceCube 
and Km3Net. If on the other hand the gamma-ray emis- 
sion is primarily leptonic, it will be accompanied by X- 
ray emission, and their extragalactic background flux in 
the X-ray band may provide an additional constraint. 
We will pursue such models and calculations in an up- 
coming publication. 

In this work we have tried, where possible, to make as- 
sumptions, which, if anything, underestimate the possi- 
ble contribution of unresolved unidentified sources to the 
EGRB. An exception to this general trend is our work- 
ing assumption that the majority of the resolved, high- 
latitude 3EG unidentified sources belong to a single class. 
It is conceivable that instead, the resolved unidentified 
sources are a collection of members of several known and 
unknown classes of gamma-ray emitters. In this case, it 
is still likely that the summed contribution of unresolved 
members of all parent classes to the diffuse background is 
significant. However, the construction of a single cumu- 
lative flux distribution from all sources and its extrapo- 
lation to lower fluxes is no longer an indicative test for 
the importance of such a contribution. 

Although we have argued that the contribution of un- 
resolved unidentified sources to the EGRB is likely to 
be important or even dominant, the cumulative emission 



spectrum derived here and presented in Fig. [3] is only an 
upper limit, as it was derived demanding that the ob- 
served EGRB is not exceeded at any energy above 100 
MeV. The contribution of unresolved unidentified sources 
is further constrained by allowing for the presence of the 
guaranteed contributions of unresolved normal galaxies 
and blazars. 

Finally, it is noteworthy that we have found no evi- 
dence for an inconsistency between the population prop- 
erties of high-latitude unidentified sources and those of 
blazars. The presence of a yet- unknown population of ex- 
tragalactic high-energy emitters among the high- latitude 
EGRET unidentified sources remains one of the most 
tantalizing possibilities in GeV astronomy and one of the 
most exciting prospects for the GLAST era. However, 
the presently available data on the spectral index dis- 
tribution and the cumulative flux distribution of these 
sources (this work), as well as their variability proper- 
ties (e.g., Nolan et al. 2003), are consistent with their 
being members of the blazar population. If indeed a 
large number of members of the blazar class are present 
among the entragalactic unidentified sources, then this 
could have potentially serious effects on our understand- 
ing of the redshift distribution of resolved blazars and 
consequently of the blazar luminosity function. 
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APPENDIX 



A. ESTIMATING THE SPECTRAL INDEX DISTRIBUTION OF EXTRAGALACTIC UNIDENTIFIED 

SOURCES 

Figure lATl shows a histogram of the spectral indices of the sources in Sample 1 (solid line). The typical measurement 
uncertainty for any single spectral index (thick solid line) is comparable with the spread of the distribution, so the 
spread of a simple binning of spectral indices might not in fact give us information about the underlying distribution 
of the spectral indices of the sources, but rather be representative of the uncertainty of each single measurement. 

Following an analysis similar to that of Venters & Pavlidou (2007) for the case of blazars, we assume that the intrinsic 
spectral index distribution (ISID) of unidentified sources can be approximated by a Gaussian, 



p(a) 



27T<7o 



exp 



(a - apf 
2al 



(Al) 



Lacking any knowledge to the contrary, this choice not only simplifies the subsequent statistical analysis, but also 
renders the spectral shape part of Eq. ([6]) analytically integrable (see Appendix). We then use a maximum- likelihood 
analysis which takes into account the individual errors of measurement of a for each source by introducing the true 
spectral indices of the sources as nuisance parameters and by marginalizing over them. The likelihood in this case, 
omitting constant normalization factors, is (Venters & Pavlidou 2007) 



C = 




exp 



A' 



1 \ " VT I 

2 ^ 



(a d - a ) 2 
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where ay and aj are the measured spectral index and the associated uncertainty of this measurement for a single 
source. Table lAll shows the maximum- likelihood ao & n d ctq for each sample. The maximum-likelihood spectral index 
distribution for sample 1 is plotted with the dashed line in Fig. IA1[ and indeed is narrower than the histogram. 

Figure IA2I shows la, 2a, and 3a contours for the likelihood of Eq. IA21 as calculated using Sample 1 data (solid 
line). The location of the maximum is indicated by x , while the corresponding maxima of the Sample 2 and Sample 3 
likelihoods are indicated by □ and A respectively. In both cases, the location of the maximum is within the la contour 
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Fig. Al. — Spectral index distribution of resolved unidentified sources in Sample 1 (all sources unidentified in 3EG and all potential 
blazar associations). Solid line: histogram of the EGRET data. Dashed line: maximum-likelihood Gaussian. Thick solid line: typical 
uncertainty of individual spectral index determination using EGRET data. 
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Fig. A2. — Likelihood contours for the intrinsic spectral index distribution parameters of Sample 1 (solid line). The maximum-likelihood 
parameters are denoted by X . The □ and A symbols denote the locations of the maximum-likelihood parameters for Sample 2 and Sample 
3 respectively. As a comparison, the likelihood contours for the sample of confidently identified blazars (as in Mattox ct al. 2001) are 
overplotted with the dashed lines. 



TABLE Al 

Maximum-likelihood parameters for the intrinsic spectral index distribution 





Sample 1 


Sample 2 


Sample 3 


ao 


2.38 


2.39 


2.40 


o-o 


0.22 


0.24 


0.19 



of Sample 1. As an additional comparison, we overplot, with the dashed lines, the 1, 2, and 3cr likelihood contours for 
the sample of the 46 confidently identified blazars of Mattox et al. (2001). The figure demonstrates that extragalactic 
unidentified sources are spectrally consistent, at the la level, with being blazars. 



10 



B. SOURCE VARIABILITY 

Certain classes of sources contributing to the gamma-ray background may exhibit significant time variability (blazars 
are a prominent example of such sources). This fact has generated some concern that the objects resolved by EGRET 
(especially the fainter ones) are preferentially objects that happened to be in their flaring state when observed, and 
hence the luminosity function as determined by the resolved objects is not an adequate description of all the unresolved 
sources contributing to the gamma-ray background. Stecker & Salamon (1996a) have addressed that uncertainty in 
the case of blazars using an explicit treatment of source variability. 

This concern is less pronounced if the fluxes used for normalizing the luminosity function of a class of sources are 
averaged over exposure times reasonably large compared to the source duty cycle. To minimize the possible impact 
of variability, we use in this study the P1234 fluxes quoted in 3EG. Sources which have been detected in individual 
viewing periods but which have no PI 234 flux (which implies that their detection significance in the cumulative map 
was not high enough) are omitted from our analysis. 

We have not made any further attempt to treat variability in this work, because the unknown parameters in the 
variability properties of sources are severe enough that any simple treatment of variability is likely to introduce at least 
as much uncertainty as it attempts to remove, while any uncertainty in the flux distribution of the source samples due 
to source variability is minor compared to the uncertainties introduced due to our ignorance concerning the nature of 
unidentified sources. We have simply used the time-averaged flux of resolved unidentified sources as representative of 
the emission properties of such sources. 



C. SPECTRAL SHAPE OF CUMULATIVE UNRESOLVED POINT SOURCE EMISSION 

The spectral shape of the cumulative emission from a population of unresolved sources with individual power- 
law photon spectra Fe oc E~ a and a Gaussian spectra index distribution p(a) is given by Eq. If p(a) is a 
Gaussian of the form of Eq. (|Aip . as assumed here, then Eq. ([6]) can be integrated analytically. Defining Iq = 
I [{kCF 8 k+1 )/(1 - k)] we can write 
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where A = a — OQ\n{E/Eo) and B Q = 2a O 0o ln(E/E ) — a^\ki(E / E^)] 2 . The usual way of plotting the isotropic 
background spec trum is E 2 Ie vs E in logarithmic scale (as in Fig. [3]). Defining X = \sx{E/Eq) and Y — \u(E 2 Ie), we 
can rewrite Eq. (|Clj) as 



Y 



£X 2 - (ao -2)X + C + ln(ao - 1 



(C2) 



where C = ln(IoE 2 ). Note that the limits of integration of Eq. [UT] should strictly extend from a = 1 (instead of — oo) 
to oo, to ensure that all F(> 100 MeV) are positive and finite. The result of Eq. (|Cip is a good approximation as long 
as ao is high enough and <7o is low enough so that p(a) has no appreciable power at a < 1. Equation (|C2|) consists of a 
parabolic convex part and a logarithmic concave part. The parabolic part dominates at low energies and is responsible 
for the characteristic convex shape of the unresolved emission from a population of sources with power-law spectra 
and a finite-spread spectral index distribution (see, e.g., Brecher & Burbidge 1972; Stecker & Salamon 1996 a,b; Pohl 
et al. 1997). The spectrum remains convex (d 2 Y/dX 2 is positive) for X < (ao — l)/o"o — !• The curvature of the 
spectrum depends on the spread of the spectral index distribution, Co (the second-order term coefficient in Eq. IC2I 
is (Tq/2). Consequently, an overestimation of the spread of the population spectral index distribution will lead to an 
overestimation of the curvature of the corresponding unresolved emission spectrum (as pointed out by Pohl et al. 1997, 
Pavlidou et al. 2006, and Venters & Pavlidou 2007). In the limit ao — > all sources have identical spectra of index 
ao, and Eq. (|C2[) yields a power-law, Y = const + (ao — 2)X or, equivalently, Ie oc E~ a °, as expected. For a < 0.5, 

Eq. (|C2[) has a local minimum, at X — (2ao — 3 — y/l — 4ctq)/(2(Tq), which corresponds to the energy at which the 
population has the smallest energy flux contribution to the (unabsorbed) isotropic background. 
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